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1.  Introduction 


Urbanization  transforms  a  natural  landscape  into  a  built  environment,  thus  affecting  the 
energy  and  water  balances  across  multiple  spatial  and  temporal  scales  [e.g.,  Georgescu  et  al., 
2009;  Lee  et  al.,  2012;  Chow  et  al,  2014].  Urban  areas  are  comprised  of  construction  materials 
(concrete,  asphalt,  gravel  cover,  etc.),  vegetation  of  different  types  (native  and  non-native,  with 
or  without  irrigation),  and  transportation  and  industrial  activities.  Each  uniquely  impacts  the 
water  and  energy  exchanges  between  the  urban  surface  and  the  overlying  atmosphere. 

Urban  water  and  energy  exchanges  can  be  particularly  important  in  desert  cities,  such  as 
Phoenix,  Arizona,  which  naturally  would  be  a  water-limited  ecosystem  of  the  Sonoran  Desert. 
However,  urban  development  leads  to  the  introduction  of  outdoor  water  use  [ Gober  et  al.,  2010; 
Volo  et  al.,  2014],  which  alters  the  water  balance.  Similarly,  the  urban  heat  island  (UHI)  effect 
alters  the  energy  balance.  Urban  water  and  energy  fluxes  between  the  surface  and  atmosphere  are 
controlled  by  land  surface  characteristics  [ Oke  et  al.,  2006], 

The  eddy  covariance  method  is  a  popular  technique  to  measure  water,  energy,  and  carbon 
fluxes  over  a  local  scale.  Recent  studies  have  used  the  eddy  covariance  method  to  measure  urban 
surface  energy  balance  [e.g.  Grimmond  and  Christen,  2012;  Chow  et  al.,  2014;  Crawford  and 
Christen,  2015].  The  eddy  covariance  footprint  is  the  surface  area  influencing  the  flux 
measurements,  and  constantly  changes  as  a  function  of  wind  direction,  atmospheric  stability,  and 
lateral  dispersion  qualities  of  the  flow  [ Schmid ,  1994].  A  major  challenge  remains  in  quantifying 
how  land  surface  characteristics  within  the  eddy  covariance  footprint  impact  the  flux 
measurements  [Lori dan  and  Grimmond,  2012] 

In  this  study,  a  mobile  eddy  covariance  platform  is  deployed  at  three  urban  settings 
located  on  Arizona  State  University  (ASU)  campuses  representing  three  urban  land  cover  types: 
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(1)  xeric  landscape  (gravel  cover  and  palo  verde  trees  with  drip-irrigation  systems  near  tall 
buildings);  (2)  high-density  urban  site  (asphalt-paved  parking  lot  near  a  high-traffic  intersection); 
and  3)  suburban  mesic  landscape  (turf  grass  with  sprinkler  irrigation  in  a  residential 
neighborhood).  The  objectives  are  to:  (a)  measure  energy,  water  and  carbon  fluxes  across 
different  land  surface  types;  (b)  quantify  the  tower  footprint  extent  to  identify  contributing 
source  areas  to  the  flux  measurements;  and  (c)  analyze  the  differential  response  to  precipitation 
and  irrigation  input  at  the  various  land  cover  sites.  This  analysis  will  improve  our  understanding 
of  how  energy  and  water  fluxes  vary  across  desert  urban  areas. 

2.  Methods 

2.1.  Study  Sites  and  Characteristics 

The  three  sites  in  this  study  are  located  on  ASU  property  (Table  1),  representing  common 
urban  landscapes  observed  in  the  Phoenix  metropolitan  area.  The  first  site  was  established  in  a 
xeric  landscape  on  the  ASU  Tempe  campus  (Fig.  la,b).  The  xeric  landscape  is  composed  of  palo 
verde  trees  with  gravel  and  bare  soil.  The  palo  verde  trees  are  regularly  irrigated  using  a  drip 
irrigation  system  and  are  approximately  3  to  5  meters  tall.  The  duration  of  the  mobile  eddy 
covariance  tower  (ECT)  deployment  at  the  palo  verde  site  (PV-ECT)  was  53  days  from  January 
20,  2015  to  March  13,  2015.  The  second  site  was  established  on  the  ASU  Tempe  campus  in  a 
parking  lot,  with  a  large  extent  of  paved  surfaces  (Fig.  lc,d).  The  duration  of  the  parking  lot 
mobile  eddy  covariance  tower  (PL-ECT)  deployment  was  43  days  from  May  19,  2015  to  June 
30,  2015.  The  third  site  was  established  near  residences  on  the  ASU  Polytechnic  campus  in 
Mesa,  Arizona,  to  represent  a  mesic  landscape  (Fig.  le,f).  The  landscape  consists  of  turf  grass, 
which  is  regularly  irrigated  using  overhead  sprinklers  three  times  per  week  and  mowed  two  times 
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per  week.  The  duration  of  the  turf  grass  mobile  eddy  covariance  tower  deployment  (TG-ECT) 
was  74  days  from  July  9,  2015  to  September  18,  2015. 

Following  Stewart  and  Oke  [2012],  each  site  is  classified  to  a  Local  Climate  Zone  (LCZ) 
to  qualitatively  describe  its  characteristics.  PV-ECT  associates  with  LCZ  C,  which  has  a  zone 
function  of  natural  scrubland  and  is  defined  as  an  open  arrangement  of  bushes,  shrubs,  and  short, 
woody  trees.  PL-ECT  is  described  with  LCZ  E,  which  is  defined  as  a  featureless  landscape  of 
rock  or  paved  cover,  few  or  no  trees  or  plants.  Its  zone  function  is  a  natural  desert  or  urban 
transportation.  TG-ECT  is  described  with  LCZ  6,  defined  as  an  open  arrangement  of  low-rise 
buildings  and  an  abundance  of  pervious  land  cover  (low  plants  and  scattered  trees). 

The  climate  of  Phoenix,  AZ  and  Mesa,  AZ  is  classified  as  a  hot  arid  subtropical  desert 
(Koppen  climate  classification  BWh).  Average  annual  temperature  is  75°F  and  71°F  in  Phoenix 
and  Mesa,  respectively.  Precipitation  is  bimodal,  occurring  predominantly  during  the  winter 
(December  to  February)  with  low  intensity  and  longer  duration  storms,  and  the  summer  (July  to 
September)  with  short,  high  intensity  storms,  which  are  representative  of  the  North  American 
monsoon  (NAM)  [Adams  and  Connie ,  1997].  Average  annual  precipitation  is  approximately  210 
mm/yr  in  Phoenix  and  approximately  235  min/yr  in  Mesa. 

2.2.  Land  Cover  Classifications 

Each  site  represents  a  different  pattern  of  distinct  urban  land  cover  features.  A  land  cover 
classification  was  perfonned  on  high-resolution  (1  foot)  orthoimagery  obtained  from  the  U.S. 
Geological  Survey.  Classifications  were  based  on  the  Red,  Green,  and  Blue  (RGB)  bands  of  the 
orthoimagery  using  a  maximum  likelihood  classification  that  utilized  training  samples  that  were 
verified  with  ground  observations.  Land  cover  was  classified  into  five  types  for  the  PV-ECT 
(Fig.  2c)  and  PL-ECT  (Fig.  3c)  sites,  and  six  types  for  the  TG-ECT  site  (Fig.  4c):  Trees;  Grass; 
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Undeveloped-Open  Land  (generally  gravel  or  bare  soil);  Pavement;  Buildings  and  Concrete;  and 
Residential  Buildings  (TG-ECT  site  only).  Buildings  and  concrete  were  classified  together  since 
both  appear  gray-white  in  the  imagery.  Trees,  grass,  and  undeveloped-open  land  represent 
pervious  surfaces,  while  pavement,  buildings  and  concrete,  and  residential  buildings  represent 
impervious  surfaces. 

2.3.  Mobile  Eddy  Covariance  Tower  Platform  and  Flux  Processing 

The  mobile  eddy  covariance  (EC)  tower  platform  consists  of  a  telescoping  tower 
mounted  on  a  trailer  that  can  extend  to  a  height  of  15  meters.  The  tower  was  set  at  different 
heights  among  the  three  different  land  surfaces  to  attempt  to  measure  fluxes  above  the  canopy  or 
building  heights  in  its  fetch,  yet  maintain  a  relatively  small  footprint  so  that  the  land  cover  of 
interest  comprised  the  majority  of  the  tower  fetch  (Table  2).  Instruments  are  mounted  to  measure 
meteorological  fluxes  between  the  urban  surface  and  atmosphere,  including:  barometric  pressure 
(CS100,  Setra);  air  temperature  and  humidity  at  three  different  heights  (near  ground  surface  ~0.3 
m,  ~0.5*height  of  EC  instruments,  and  -height  of  EC  instruments),  (HMP155A,  Vaisala); 
incoming  solar  radiation  at  5  m  height  (CMP3,  Kipp  &  Zonen);  net  incoming  and  outgoing 
longwave  and  shortwave  radiation  at  5  m  height  (CNR4,  Kipp  &  Zonen);  incoming  and  outgoing 
photosynthetic  active  radiation  (PAR)  (SQ-1 10,  Apogee  Instruments  Inc.);  incoming  and 
outgoing  shortwave  radiation  (SP-1 10,  Apogee  Instruments  Inc.);  C02  and  H20  concentrations 
using  an  open-path  Infrared  Gas  Analyzer  (LI7500A,  LI-COR);  and  three-dimensional  wind 
velocity  (CSAT3,  Campbell  Sci.)  at  7,  8  or  9  m  height  aligned  to  the  dominant  wind  direction  of 
the  season  as  determined  by  wind  rose  diagrams  from  nearby  meteorological  stations  (Table  2). 
Soil  conditions  are  evaluated  with  measurements  of:  surface  temperature  (SI-1 1 1,  Apogee 
Instruments  Inc.);  soil  moisture  at  5,  15,  and  50  cm  depths  (CS616,  Campbell  Sci.);  soil 
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temperature  at  2  and  4  cm  depths  (TCAV,  Campbell  Sci.);  and  soil  heat  flux  at  5  cm  depth 
(HFP01-SC,  Hukseflux).  Soil  moisture,  soil  temperature  and  soil  heat  flux  were  not  measured  at 
the  PL-ECT  site  as  the  tower  was  measuring  a  paved  surface.  EC  data  were  sampled  at  a  20  Hz 
frequency  for  the  PV-ECT  site  and  at  10  Hz  frequency  for  the  PL-ECT  and  TG-ECT  sites  and 
recorded  by  a  datalogger  (CR5000,  Campbell  Sci.  or  CR3000,  Campbell  Sci.). 

Fluxes  were  processed  using  EdiRE  (University  of  Edinburgh),  which  includes  correction 
for  fluctuations  in  stability  [ Foken  et  al.,  2006]  and  density  [  Webb  et  al.,  1980],  sonic 
temperature  use  to  calculate  sensible  heat  flux,  rotating  the  coordinate  frame  to  set  the  mean 
vertical  wind  speed  to  zero  [  Wilczak  et  al.,  2001],  and  removing  signal  lag  in  gas  concentrations 
[Massman,  2001].  Flux  data  was  fdtered  to  exclude  time  periods  when  there  was  precipitation, 
the  wind  direction  was  180°  ±10°  from  the  direction  at  which  the  EC  instruments  were  mounted 
due  to  possible  interference  from  the  tower  setup,  and  for  outliers  greater  than  3  standard 
deviations.  Other  measurements  were  recorded  by  a  datalogger  (CR5000,  Campbell  Sci.  or 
CR3000,  Campbell  Sci.)  as  averages  every  30  minutes. 

2.4.  Footprint  Analysis 

Eddy  covariance  measurements  were  used  to  estimate  the  footprint,  or  source  area,  of  the 
measured  turbulent  fluxes  [Gash,  1985;  Schmid,  1997].  The  analytical  model  from  Kormann  and 
Meixner  [2001]  was  used  to  calculate  the  footprint  for  an  area  of  600  meters  by  600  meters, 
centered  at  each  ECT.  The  horizontal  resolution  of  the  footprint  contributions  was  set  to  3 
meters,  which  is  less  than  the  measurement  heights  [  Van  de  Boer  et  al.,  2013;  Anderson  and 
Vivoni,  2015].  The  footprint  calculations  were  processed  using  EDiRE,  and  depend  on 
measurement  height,  mean  horizontal  wind  velocity,  wind  direction,  friction  velocity,  Monin- 
Ohbukov  stability,  and  the  standard  deviation  of  the  lateral  wind  velocity.  Lower  limits  of  wind 
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velocity,  friction  velocity,  stability  and  sensible  heat  flux  were  set  to  ensure  footprints  were 
calculated  during  turbulent  daytime  conditions.  Footprints  were  estimated  at  30  minute  time 
intervals  and  then  averaged  at  a  daily  scale.  The  daily  averaged  footprints  were  then  used  to 
derive  a  unique  deployment  period  averaged  footprint  for  each  ECT. 

Source  area  represents  the  area  around  each  ECT  that  contributes  to  the  flux 
measurements.  In  this  study,  the  50%  source  area  and  the  80%  source  area  are  evaluated.  The 
50%  source  area  is  defined  as  the  area  that  is  responsible  for  50%  of  the  surface  influence,  with 
regards  to  the  flux  measurements  [ Schmid ,  1994].  The  area  that  accounts  for  80%  of  the  surface 
influence  is  referred  to  as  the  80%  source  area. 

3.  Results  and  Discussion 

3.1.  Footprint  Representation  of  Land  Cover 

Land  cover  classification  was  performed  using  1  foot  resolution  orthoimagery  centered  at 
each  ECT.  To  quantify  the  land  cover  types  that  were  contributing  to  the  flux  measurements,  the 
50%  source  areas  were  calculated  for  each  ECT.  The  80%  source  areas  were  calculated  for  PV- 
ECT  and  PL-ECT,  however  the  600  meter  by  600  meter  box  centered  at  the  TG-ECT  did  not 
account  for  at  least  80%  of  the  flux  measurements  at  the  site.  This  could  be  improved  in  a  future 
effort  by  defining  a  larger  fetch.  The  percent  of  area  that  was  accounted  for  in  the  300  meter 
fetch  (or  600  meter  by  600  meter  box)  is  96.2%,  90.1%  and  61.1%  for  the  PV-ECT,  PL-ECT, 
and  TG-ECT,  respectively.  As  a  result,  the  50%  source  areas  are  used  to  quantitatively  compare 
land  class  distribution  among  the  three  sites  in  this  study  (Table  3). 

The  PV-ECT  site  (Fig.  2b, c)  has  a  relatively  small  50%  source  area  (1,917  m2).  Trees, 
specifically  palo  verde  trees,  dominate  the  land  cover  at  56%.  The  next  dominant  land  cover  is 
undeveloped  (bare  soil  or  gravel)  at  37%.  Pavement,  buildings  and  concrete,  the  impervious 
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materials,  account  for  7.8%  together.  When  a  larger  source  area  is  considered  (80%  source  area), 
the  land  cover  classes  are  more  evenly  distributed,  with  the  footprint  extending  beyond  the  palo 
verde  patch  to  surrounding  buildings,  concrete,  and  pavement.  The  PL-ECT  site  (Fig.  3b, c)  has  a 
larger  50%  source  area  (9,243  m  ).  The  most  dominant  land  class  is  pavement,  as  expected,  at 
83%.  Undeveloped  regions  account  for  10%.  Trees  and  grass  make  up  4.5%  with  buildings  and 
concrete  at  2.5%.  Again,  the  land  cover  classes  are  more  evenly  distributed  when  a  larger  source 
area  is  considered,  however  pavement  still  dominates  at  56%  in  the  80%  source  area.  The  TG- 
ECT  site  (Fig.  4b, c)  is  the  most  evenly  distributed.  It  has  the  largest  50%  source  area  of  the  three 
sites  at  14,346  m  .  The  land  cover  classifications  are  much  more  evenly  distributed.  The 
dominant  land  class  is  undeveloped  (bare  or  gravel)  land  at  36.5%  followed  by  grass  cover  at 
23.6%.  Trees  account  for  18%,  pavement  is  1 1.5%  and  then  residential  buildings  comprise  8%. 
Building  and  concrete  have  the  smallest  areal  percentage  at  2.6%.  Pervious  land  cover  types 
comprise  approximately  78%  of  the  50%  source  area. 

3.2.  Hydrologic  Dynamics  Across  Variable  Seasons 

The  PV-ECT  deployment  period  was  in  the  winter  time,  however  it  had  the  greatest 
amount  of  total  rainfall  (Table  4).  The  largest  storm  also  had  the  longest  duration  (~1  day)  and 
accumulated  23.8  mm  between  18:00  on  1/29/2015  to  20:00  on  1/30/2015.  In  general  the  PV- 
ECT  measured  lower  temperatures  (surface  and  air)  due  to  the  season  (Fig.  5a).  Surface 
temperature  had  the  highest  variation  and  had  much  higher  maximum  daily  temperature  values. 
The  surface  temperature  was  measuring  the  temperature  of  the  inter-canopy  spaces  which 
consists  of  gravel  and  bare  soil.  Incoming  solar  and  net  radiation  were  also  the  lowest  of  the 
three  deployments  because  of  the  season  (Fig.  5b).  During  rainfall  events,  incoming  solar  and  net 
radiation  decreased  due  to  cloud  cover.  The  PV-ECT  generally  also  had  high  sensible  heat  fluxes 
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(Fig.  5c).  There  was  one  large  increase  in  latent  heat  flux  after  the  day  long  storm  event  at  the 
end  of  January.  This  is  expected  as  the  moisture  input  was  large  and  palo  verde  trees  represent  a 
large  areal  fraction  of  the  footprint.  A  slight  increase  in  latent  heat  flux  was  also  observed  after 
an  early  March  rainfall  event.  Soil  moisture  was  measured  at  two  depths,  5  cm  and  50  cm  (Fig. 
5d).  The  shallow  depth  has  a  slight  response  every  day,  which  is  assumed  to  be  related  to  the  drip 
irrigation,  and  a  large  response  to  precipitation  input.  The  50  cm  depth  remained  relatively 
constant  during  the  deployment  period.  This  is  expected  because  during  the  installation  of  the 
sensor,  a  black  mesh  isolating  the  lower  soil  was  discovered  at  -20  cm.  The  mesh  will  limit  the 
amount  of  water  that  can  infiltrate.  This  may  also  explain  the  increased  latent  heat  flux  after  the 
long  storm  even  in  January.  A  large  precipitation  input  coupled  with  the  limited  ability  to 
infiltrate  water  could  result  in  high  surface  evaporation  rates. 

The  PL-ECT  deployment  period  occurred  from  late  May  through  June,  which  is  typically 
the  driest  and  hottest  time  of  the  year.  Surface  temperature  was  very  high  (Fig.  6a).  The  sensor 
was  measuring  the  gravel  that  mobile  ECT  platform  sat  on,  however  temperatures  measured 
were  still  very  high,  with  peaks  up  to  68°C.  Air  temperature  was  relatively  high,  due  to  the  time 
of  year.  The  air  temperature  measurement  at  the  ‘high’  height,  approximately  9  m,  had  the  lowest 
readings,  with  the  air  temperature  measurement  closest  to  the  surface  giving  the  highest  readings. 
Incoming  solar  radiation  and  net  radiation  were  relatively  high  (Fig.  6b).  During  the  precipitation 
event  in  early  June,  temperatures,  incoming  solar  and  net  radiation  all  decreased.  Latent  heat  flux 
also  increased  slightly  following  the  early  June  rainfall  event  (Fig.  6c).  However,  sensible  heat 
flux  was  dominant  for  the  entire  deployment  period,  largely  due  to  concrete  representing  a  large 
majority  of  the  ECT  footprint.  Ground  heat  flux  could  not  be  measured  at  the  site.  Interestingly, 
there  is  not  a  large  change  in  any  of  the  measured  variables  to  the  precipitation  events  that 
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occurred  on  June  27  and  June  29.  This  may  be  because  both  events  occurred  at  night,  and  both 
were  short  with  high  intensity.  It  is  possible  that  summer-type  storms  have  a  small  temporal 
impact  on  pavement  land  cover  with  regards  to  flux  measurements. 

The  TG-ECT  deployment  occurred  during  the  NAM  season,  from  mid-July  to  mid- 
September.  Surprisingly  however,  there  was  relatively  little  rainfall  recorded  during  the  period 
(Table  4).  In  contrast  to  the  other  two  sites,  surface  temperature  was  lower  than  air  temperature, 
corresponding  to  the  irrigated  turf  grass  (Fig.  7a).  As  the  air  temperature  was  taken  at  increasing 
heights,  the  air  temperature  measurements  also  increased.  Incoming  solar  and  net  radiation  were 
relatively  high,  as  expected  for  the  season  (Fig.  7b).  Also  in  contrast  to  the  other  two  sites,  latent 
heat  flux  was  dominant  (Fig.  7c).  This  is  expected  as  the  grass  is  regularly  irrigated  (three  times 
per  week).  Latent  heat  flux  appears  to  slightly  increase  after  rainfall  events.  Soil  moisture  was 
measured  at  5  cm,  15  cm,  and  50  cm  depths  (Fig.  7d).  All  of  the  soil  moisture  values  were  high. 
The  5  cm  and  1 5  cm  soil  moisture  appears  to  vary  ~3  times  per  week,  mimicking  the  irrigation 
schedule.  A  deeper  moisture  pulse  was  observed  with  the  rainfall  event  on  8/1 1/2015,  which  had 
1.7  mm  of  measured  precipitation.  Otherwise,  soil  moisture  at  50  cm  was  fairly  constant. 

The  diurnal  cycle  for  energy  flux  components  was  averaged  at  30  minute  time  intervals 
for  each  site  over  each  site’s  deployment  time  period  (Fig.  9).  PL-ECT  has  the  highest  net 
radiation,  which  is  expected  given  the  time  of  year  for  the  deployment.  At  nighttime,  PL-ECT 
also  has  the  most  negative  net  radiation.  This  is  likely  due  to  the  outgoing  longwave  radiation 
correlated  with  the  high  surface  temperature  of  the  pavement.  Sensible  heat  flux  is  fairly  similar 
between  the  PV-ECT  and  the  PL-ECT.  There  is  a  slight  lag  of  peak  sensible  heat  flux  at  PL- 
ECT,  which  occurs  around  14:00.  Sensible  heat  flux  is  lowest  at  TG-ECT.  TG-ECT  has  the 
highest  latent  heat  flux,  which  peaks  around  12:00.  PL-ECT  has  the  lowest  latent  heat  flux, 
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which  peaks  slightly  around  13:00.  Ground  heat  flux  patterns  are  similar  between  PV-ECT  and 
TG-ECT.  Average  rainfall  is  computed  as  the  average  amount  of  rainfall  that  fell  during  that  30 
minute  time  period  over  the  duration  of  the  deployment.  The  one  day  long  precipitation  event 
during  the  PV-ECT  deployment  is  evident  in  this  calculation  (Fig.  9a)  because  there  is  a 
contribution  for  almost  every  30  minute  period. 

To  evaluate  the  energy  balance,  the  30  minute  averaged  flux  data  was  plotted  with  net 
radiation  minus  ground  heat  flux  on  the  x-axis,  and  sensible  heat  flux  plus  latent  heat  flux  on  the 
v-axis.  A  linear  regression  was  fit  to  each  plot  and  the  y-intercept  and  slope  were  calculated, 
following  this  equation:  LE  +H=m(Rn-G)+b.  For  PV-ECT,  v  =  0.53x  +  27  with  R2  =  0.90.  The 

result  of  the  linear  regression  to  evaluate  the  energy  balance  lit  for  TG-ECT  is:  y  =  0.72x  +  33 
with  R  =  0.89.  This  indicates  that  TG-ECT  site  has  a  better  energy  balance  closure. 

3.3.  Comparison  of  Urban  Surface  Fluxes 

To  better  compare  the  sites  to  one  another,  different  energy  flux  ratios  were  evaluated. 
Evaporative  fraction  (EF)  is  computed  at  noon  time  of  each  day  as  EF  =  LE/(LE+H),  where  LE 
is  latent  heat  flux  and  El  is  sensible  heat  flux.  Bowen  ratio  ( B )  is  computed  at  noon  time  of  each 
day  as  B  =  H/LE.  An  average  was  taken  of  the  computed  EF  and  B  over  each  site’s  deployment 
period  and  they  are  summarized  in  Table  5.  Ratios  of  sensible  heat  to  net  radiation  minus  ground 
heat  flux  and  latent  heat  to  net  radiation  minus  ground  heat  flux  were  also  computed  at  a  daily 
average  time  step  (Fig.  10),  and  will  be  referred  to  as  the  sensible  heat  ratio  and  latent  heat  ratio, 
respectively.  The  average  over  all  of  the  days  in  each  site’s  deployment  period  was  calculated 
(Table  5).  On  average,  evaporative  fraction  was  highest  at  TG-ECT,  followed  by  PV-ECT  and 
then  PL-ECT.  This  is  expected  due  to  the  high  measured  latent  heat  fluxes  at  the  turf  grass  site, 
and  the  very  low  latent  heat  fluxes  at  the  parking  lot  site.  Bowen  ratio  is  ~3  at  PV-ECT  and  ~8  at 
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PL-ECT,  which  are  representative  of  a  semi-arid  to  desert  environment  because  sensible  heat 
fluxes  are  high.  Average  B  is  0.53  at  TG-ECT,  which  means  that  latent  heat  flux  is  dominant. 
When  evaluating  the  time  series  of  the  flux  ratios,  the  rainfall  events  that  occurred  at  PV-ECT 
are  apparent.  Sensible  heat  ratio  is  generally  higher  except  for  the  few  days  around  each  rainfall 
event.  At  PL-ECT,  both  ratios  are  likely  underestimated  due  to  no  ground  heat  flux 
measurement,  therefore  G  was  assumed  to  equal  zero  for  the  ratio  calculation.  The  sensible  heat 
ratio  is  always  higher  than  the  latent  heat  ratio.  The  difference  between  the  two  is  fairly  constant, 
but  decreases  during  the  rainfall  event  in  early  June.  TG-ECT  observes  an  opposite  trend,  as 
latent  heat  ratio  is  generally  always  higher  than  sensible  heat  ratio.  However,  the  difference 
between  the  two  is  much  more  variable. 

To  further  evaluate  how  the  surfaces  responded  to  precipitation  input,  EF,  B ,  sensible 
heat  and  latent  heat  ratios  were  calculated  for  ‘Wet’,  ‘Dry’,  ‘Wet  Daytime’  (abbreviated  as  ‘Wet- 
Day’),  and  ‘Dry  Daytime’  (abbreviated  as  ‘Dry-Day’)  at  averaged  day  periods  (Table  6).  ‘Wet’ 
days  are  defined  as  any  day  where  precipitation  occurred  or  a  day  where  precipitation  occurred 
within  the  last  3  days.  Essentially  averages  were  computed  for  the  day  of  rainfall  plus  the  three 
subsequent  days  after  rainfall.  ‘Dry’  days  are  defined  as  all  other  days  during  the  deployment 
period  that  did  not  meet  that  criteria.  ‘Wet-Day’  time  steps  are  defined  as  days  where 
precipitation  occurs  during  day-time  hours  and  three  days  after  such  a  precipitation  event.  ‘Dry- 
Day’  is  defined  as  all  other  days  during  the  deployment  period  that  do  not  meet  the  ‘Wet-Day’ 
criteria.  In  summary,  PV-ECT  has  18  and  33  wet  and  dry  days,  with  13  and  38  wet-day  and  dry- 
day  time  periods.  PL-ECT  has  9  and  32  wet  and  dry  days,  with  6  and  35  wet-day  and  dry-day 
time  periods,  and  TG-ECT  has  17  and  53  wet  and  dry  days,  with  7  and  63  wet-day  and  dry-day 
time  periods.  By  differentiating  wet  and  dry  days,  the  patterns  observed  are  strengthened.  All 
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sites  have  an  increased  EF  on  wet  days  as  compared  to  dry  days,  with  the  difference  increasing 
the  most  at  PL-ECT.  The  difference  between  EF  is  larger  at  all  3  sites  when  evaluating  only  day¬ 
time  rainfall  events.  There  is  not  a  significant  difference  in  computed  B  at  PV-ECT  when 
evaluating  all  rainfall  events,  however  there  is  when  evaluating  day-time  rainfall  events,  with  a 
lower  B  for  wet-day  day  periods.  TG-ECT  actually  averages  a  higher  B  for  wet  periods  when  all 
rainfall  events  are  considered,  but  only  day-time  rainfall  events  lead  to  a  lower  B  for  wet-days. 
Sensible  heat  ratios  are  lower  for  wet  days  and  higher  for  dry  days  at  PV-ECT  and  TG-ECT. 
Latent  heat  ratios  are  higher  or  wet  days  and  lower  for  dry  days  at  PV-ECT.  The  ratios  appear  to 
be  largely  dependent  on  moisture  availability  at  PV-ECT.  Interestingly,  the  sensible  heat  and 
latent  heat  ratios  are  both  higher  for  wet  days  and  wet-day  daily  periods  compared  to  the  dry 
values  at  PL-ECT.  It  is  important  to  note  that  G  is  not  measured  at  the  PL-ECT,  so  the 
calculation  is  a  fraction  of  sensible  heat  and  latent  heat  to  net  radiation.  TG-ECT  is  insensitive  to 
latent  heat  flux,  regardless  of  the  days  considered,  which  is  likely  due  to  the  constant  irrigation 
input.  In  general,  wet  time  periods  result  in  higher  evaporative  fractions,  lower  Bowen  ratios, 
lower  sensible  heat  ratios,  and  higher  latent  heat  ratios.  These  trends  are  more  pronounced  when 
only  day-time  precipitation  events  are  considered,  suggesting  that  the  timing  of  precipitation  has 
an  influence  on  how  the  surfaces  respond  to  water  and  energy  fluxes. 

Carbon  dioxide  fluxes  were  calculated  at  30  minute  averages,  and  then  averaged  over  a 
daily  time  step  for  each  ECT  site  to  observe  general  magnitudes  (Fig.  1 1).  Both  PV-ECT  and  PL- 
ECT  have  positive  CO?  fluxes,  indicating  that  CO2  is  being  released  into  the  atmosphere.  This  is 
likely  a  result  of  urban  CO2  sources.  The  TG-ECT  site  measures  CO2  fluxes  that  are  generally 
negative,  indicating  the  land  surface  may  sequester  carbon. 
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4.  Summary  and  Conclusions 

Urban  water  and  energy  exchanges  between  the  surface  and  atmosphere  are  controlled  by 
land  surface  characteristics.  In  this  study,  water  and  energy  fluxes  were  evaluated  at  three 
different  urban  land  cover  representations.  Comparisons  between  the  three  sites  reveal: 

(1)  Surface  temperature  is  higher  than  air  temperature  for  PV-ECT  and  PL-ECT, 
whereas  surface  temperature  is  less  than  air  temperature  at  TG-ECT. 

(2)  Sensible  heat  flux  dominates  at  PV-ECT  and  PL-ECT,  while  latent  heat  flux 
dominates  at  TG-ECT.  This  is  likely  due  to  above  ground  sprinkler  irrigated  grass 
comprising  a  sizeable  fraction  of  the  TG-ECT  footprint. 

(3)  PV-ECT  site  responds  to  long  precipitation  events  with  a  spike  in  latent  heat  flux. 
The  PL-ECT  site  appears  to  be  relatively  insensitive  to  rainfall  events,  however 
the  sample  size  of  precipitation  events  is  small. 

(4)  Precipitation  influences  B,  EF,  sensible  and  latent  heat  ratios  over  all  three  land 
surfaces.  Changes  are  more  noticeable  during  day-time  precipitation  events, 
which  may  suggest  that  precipitation  timing  influence  the  surface  response  to 
energy  and  water  fluxes. 

The  three  different  site  deployments  occurred  during  different  times  of  the  year,  which 
limits  the  ability  to  compare  fluxes  from  one  site  to  another.  Further  data  analysis  can  be  useful 
to  better  understand  difference  between  the  three  sites,  including  a  direct  comparison  of  the  three 
sites  with  continuous  data  from  an  urban  tower  in  Phoenix,  as  described  by  Chow  et  al.  [2014], 
and  a  natural  landscape  tower  in  the  Sonoran  Desert,  as  described  by  Pierini  et  al.  [2014],  As 
desert  urban  areas  continue  to  grow,  it  is  important  to  improve  our  understanding  of  how  the 
energy  and  water  balances  vary  across  urban  land  cover. 
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Figure  Captions 

Figure  1.  Three  study  site  locations  including  (a)  palo  verde  eddy  covariance  tower  (PV-ECT)  at 
ASU  Tempe  campus  in  a  (b)  xeric  landscape  consisting  of  palo  verde  trees  with  gravel  surface 
that  is  drip  irrigated,  (c)  parking  lot  eddy  covariance  tower  (PL-ECT)  at  ASU  Tempe  campus  on 
a  (d)  gravel  area  within  the  paved  parking  lot,  near  a  high  traffic  intersection,  and  (e)  turf  grass 
eddy  covariance  tower  (TG-ECT)  near  residential  housing  at  ASU  Polytechnic  campus  in  Mesa, 
AZ  in  a  (f)  mesic  landscape  consisting  of  regularly  irrigated  turf  grass.  The  0.30  m  aerial 
photographs  in  (b)  are  from  USGS  high  resolution  orthoimagery. 

Figure  2.  (a)  Ortho  image  centered  at  PV-ECT,  with  (b)  the  flux  footprint  of  the  ECT, 
representing  50%  and  80%  source  area  with  (c)  land  cover  classification  consisting  of:  trees, 
grass,  undeveloped  -  open  land  (gravel  or  bare  soil),  pavement,  and  buildings  and  concrete. 

Figure  3.  (a)  Orthoimage  centered  at  PL-ECT,  with  (b)  the  flux  footprint  of  the  ECT, 
representing  50%  and  80%  source  area  with  (c)  land  cover  classification  consisting  of:  trees, 
grass,  undeveloped  -  open  land  (gravel  or  bare  soil),  pavement,  and  buildings  and  concrete. 

Figure  4.  (a)  Ortho  image  centered  at  TG-ECT,  with  (b)  the  flux  footprint  of  the  ECT, 
representing  50%  source  area  with  (c)  land  cover  classification  consisting  of:  trees,  grass, 
undeveloped  -  open  land  (gravel  or  bare  soil),  pavement,  buildings  and  concrete,  and  residential 
buildings. 

Figure  5.  Meteorological  variables  measured  at  PV-ECT,  including  (a)  air  temperature  at  three 
different  heights  (low,  medium,  and  high)  and  surface  temperature,  (b)  incoming  solar  radiation 
(Is)  and  net  radiation  (Rn),  (c)  sensible  (H),  latent  (LE)  and  ground  (G)  heat  flux  and  (d) 
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volumetric  soil  moisture  at  5  cm  and  50  cm  depths.  Rainfall  is  plotted  as  the  total  accumulation 
[mm]  per  30  minute  time  period. 

Figure  6.  Meteorological  variables  measured  at  PL-ECT,  including  (a)  air  temperature  at  three 
different  heights  (low,  medium,  and  high)  and  surface  temperature,  (b)  incoming  solar  radiation 
(Is)  and  net  radiation  (Rn),  and  (c)  sensible  (H)  and  latent  (LE)  heat  flux.  Rainfall  is  plotted  as 
the  total  accumulation  [mm]  per  30  minute  time  period. 

Figure  7.  Meteorological  variables  measured  at  TG-ECT,  including  (a)  air  temperature  at  three 
different  heights  (low,  medium,  and  high)  and  surface  temperature,  (b)  incoming  solar  radiation 
(Is)  and  net  radiation  (Rn),  (c)  sensible  (H),  latent  (LE)  and  ground  (G)  heat  flux  and  (d) 
volumetric  soil  moisture  at  5  cm,  15  cm  and  50  cm  depths..  Rainfall  is  plotted  as  the  total 
accumulation  [mm]  per  30  minute  time  period. 

Figure  8.  Average  diurnal  cycle  for  energy  balance  fluxes  [net  radiation  (Rn),  sensible  heat  flux 
(H),  latent  heat  flux  (LE),  and  ground  heat  flux  (G)]  computed  at  30  minute  time  steps  over  each 
ECT  deployment  duration  at  (a)  PV-ECT  (b)  PL-ECT  and  (c)  TG-ECT.  Average  rainfall  is 
plotted  and  computed  as  the  average  rainfall  amount  for  each  30  minute  period  during  a  day. 

Figure  9.  Sensible  and  latent  heat  ratios  to  net  radiation  minus  ground  heat  radiation,  (H/(Rn-G) 
and  LE/(Rn-G),  respectively),  computed  as  an  average  daily  value  for  (a)  PV-ECT  (b)  PL-ECT 
and  (c)  TG-ECT. 

Figure  10.  Daily  averaged  carbon  dioxide  flux  values  at  (a)  PV-ECT  (b)  PL-ECT  and  (c)  TG- 
ECT. 
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Table  Captions 

Table  1.  General  characteristics  for  the  three  study  sites,  including  land  cover  description, 
location,  and  elevation. 

Table  2.  Eddy  covariance  deployment  specifications  for  the  three  study  sites,  including  angle 
and  height  the  high  frequency  sensors  were  mounted  at,  speed  at  which  EC  measurements  were 
being  taken,  and  each  deployment  duration. 

Table  3.  Land  cover  percentages  detennined  within  50%  and  80%  source  areas  for  each  site. 
Area  values  in  square  meters  for  50%  and  80%  source  area  at  each  site  and  the  percent  of  the 
flux  that  falls  within  the  specified  300  m  fetch. 

Table  4.  Dates  of  rainfall  events  that  occurred  during  the  deployment  period  at  each  tower,  with 
approximate  start  and  ending  times  and  total  accumulation. 

Table  5.  Average  evaporative  fraction,  Bowen  ratio,  sensible  heat  flux  fraction  of  net  radiation 
minus  ground  heat  flux,  and  latent  heat  flux  fraction  of  net  radiation  minus  ground  heat  flux  for 
the  three  sites.  Evaporative  fraction  and  Bowen  ratio  were  calculated  at  noon  each  day,  while  the 
daily  average  flux  values  were  used  to  find  the  fraction  values.  Averages  were  then  computed 
over  each  ECT  deployment  period. 

aPL-ECT  had  no  G  measurement,  therefore  G  was  assumed  to  equal  zero. 

Table  6.  Evaporative  fraction,  Bowen  ratio,  sensible  heat  flux  fraction  of  net  radiation  minus 
ground  heat  flux,  and  latent  heat  flux  fraction  of  net  radiation  minus  ground  heat  flux  for  the 
three  sites.  Evaporative  fraction  and  Bowen  ratio  were  calculated  at  noon  each  day,  while  the 
daily  average  flux  values  were  used  to  find  the  fraction  values.  Averages  were  then  computed 
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over  wet  days,  defined  as  a  day  where  rainfall  occurred  plus  three  days  after,  and  dry  days, 
defined  as  all  days  that  are  not  classified  as  wet.  Averages  were  also  computed  for  rainfall  events 
that  occurred  during  the  day  time  (Wet-Day)  including  three  days  after,  where  Dry-Day  is  all 
other  days  that  were  not  classified  as  Wet-Day. 

aPL-ECT  had  no  G  measurement,  therefore  G  was  assumed  to  equal  zero. 
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Table  1.  General  characteristics  for  the  three  study  sites,  including  land  cover  description, 
location,  and  elevation. 


Site 

Land  Cover 

UTM 

Easting 

(m) 

UTM 

Northing  (m) 

Latitude 

Longitude 

Elevation 

(m) 

PV-ECT 

Palo  Verde 

413797 

3698213 

33.4198 

-111.9272 

354 

PL-ECT 

Pavement 

412725 

3698373 

33.4212 

-111.9387 

356 

TG-ECT 

Turf  Grass 

436646 

3686041 

33.3116 

-111.6806 

411 

Table  1 
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Table  2.  Eddy  covariance  deployment  specifications  for  the  three  study  sites,  including  angle 
and  height  the  high  frequency  sensors  were  mounted  at,  speed  at  which  EC  measurements  were 
being  taken,  and  each  deployment  duration. 


Site 

Angle 

(deg) 

Measurement 
Height  (m) 

Speed 

Start  Time 

End  Time 

Days 

PV-ECT 

21 

7 

20  Hz 

1/20/2015  12:00 

3/13/2015  8:30 

53 

PL-ECT 

227 

9 

10  Hz 

5/19/2015  15:00 

6/30/2015  6:00 

43 

TG-ECT 

230 

8 

10  Hz 

7/9/2015  13:00 

9/18/2015  8:30 

74 

Table  2 
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Table  3.  Land  cover  percentages  detennined  within  50%  and  80%  source  areas  for  each  site. 
Area  values  in  square  meters  for  50%  and  80%  source  area  at  each  site  and  the  percent  of  the 
flux  that  falls  within  the  specified  300  m  fetch. 


PV-ECT 

PL-1 

ECT 

TG-ECT 

50%  Source 
Area 

80%  Source 
Area 

50%  Source 
Area 

80%  Source 
Area 

50%  Source 
Area 

Trees 

55.6% 

37.0% 

4.3% 

5.5% 

17.9% 

Grass 

0.0% 

0.2% 

0.2% 

0.7% 

23.6% 

Undeveloped 

36.6% 

27.9% 

10.0% 

14.8% 

36.5% 

Pavement 

5.6% 

26.5% 

82.9% 

56.0% 

11.5% 

Buildings/Concrete 

2.2% 

26.5% 

2.5% 

23.0% 

2.6% 

Residential 

0.0% 

0.0% 

0.0% 

0.0% 

7.9% 

Area  [m2] 

1,917 

15,066 

9,243 

76,158 

%  in  300  m  fetch 

96.2% 

90.1% 

Table  3 
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Table  4.  Dates  of  rainfall  events  that  occurred  during  the  deployment  period  at  each  tower,  with 
approximate  start  and  ending  times  and  total  accumulation. 


PV-ECT 

PL-ECT 

TG-ECT 

Date 

Time 

P  [mm] 

Date 

Time 

P  [mm] 

Date 

Time 

P  [mm] 

1/26/2015 

20:30-23:30 

0.6 

6/4/2015 

14:00-14:30 

0.3 

7/18/2015 

19:30-22:30 

1.5 

1/27/2015 

9:30 

0.1 

6/5/2015 

5:30-10:00 

2.5 

7/22/2015 

9:00 

0.3 

1/29/2015 

18:00-23:30 

5.3 

6/5/2015 

16:00-17:30 

0.7 

8/11/2015 

10:00-10:30 

0.4 

1/30/2015 

0:00-20:00 

18.5 

6/5/2015 

22:00-23:00 

1 

8/11/2015 

14:30 

0.1 

1/31/2015 

11:00-12:00 

1.2 

6/6/2015 

4:30-6:00 

1 

8/11/2015 

22:00-22:30 

1.6 

2/23/2015 

19:30-21:00 

0.7 

6/27/2015 

19:30-20:00 

5.7 

8/27/2015 

17:30-18:00 

1.1 

3/2/2015 

8:30-16:00 

11.9 

6/29/2015 

23:00 

4 

8/29/2015 

22:00-22:30 

0.3 

3/3/2015 

14:00 

0.3 

Table  4 
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Table  5.  Average  evaporative  fraction,  Bowen  ratio,  sensible  heat  flux  fraction  of  net  radiation 
minus  ground  heat  flux,  and  latent  heat  flux  fraction  of  net  radiation  minus  ground  heat  flux  for 
the  three  sites.  Evaporative  fraction  and  Bowen  ratio  were  calculated  at  noon  each  day,  while  the 
daily  average  flux  values  were  used  to  find  the  fraction  values.  Averages  were  then  computed 
over  each  ECT  deployment  period. 

aPL-ECT  had  no  G  measurement,  therefore  G  was  assumed  to  equal  zero. 


Site 

Evaporative  Fraction  H 

Bowen  Ratio  [-] 

H/(Rn-G)  [-1 

LE/(Rn-G)  [-1 

PV-ECT 

0.27 

3.05 

0.54 

0.36 

PL-ECT3 

0.16 

8.01 

0.55 

0.20 

TG-ECT 

0.61 

0.53 

0.30 

0.67 

Table  5 
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Table  6.  Evaporative  fraction,  Bowen  ratio,  sensible  heat  flux  fraction  of  net  radiation  minus 
ground  heat  flux,  and  latent  heat  flux  fraction  of  net  radiation  minus  ground  heat  flux  for  the 
three  sites.  Evaporative  fraction  and  Bowen  ratio  were  calculated  at  noon  each  day,  while  the 
daily  average  flux  values  were  used  to  find  the  fraction  values.  Averages  were  then  computed 
over  wet  days,  defined  as  a  day  where  rainfall  occurred  plus  three  days  after,  and  dry  days, 
defined  as  all  days  that  are  not  classified  as  wet.  Averages  were  also  computed  for  rainfall  events 
that  occurred  during  the  day  time  (Wet-Day)  including  three  days  after,  where  Dry-Day  is  all 
other  days  that  were  not  classified  as  Wet-Day. 

aPL-ECT  had  no  G  measurement,  therefore  G  was  assumed  to  equal  zero. 


Evaporative  Fraction  [-] 

Bowen  Ratio  [-] 

Condition 

Wet 

Dry 

Wet-Day 

Dry-Day 

Wet 

Dry 

Wet-Day 

Dry-Day 

PV-ECT 

0.29 

0.26 

0.32 

0.25 

3.0 

3.1 

2.0 

3.3 

PL-ECT 

0.21 

0.15 

0.22 

0.15 

4.1 

9.1 

3.6 

8.8 

TG-ECT 

0.62 

0.61 

0.69 

0.61 

0.65 

0.49 

0.47 

0.53 

H/(Rn-G)  [-] 

LE/(Rn-G)  [-] 

1 

Condition 

Wet 

Dry 

Wet-Day 

Dry-Day 

Wet 

Dry 

Wet-Day 

Dry-Day 

PV-ECT 

0.45 

0.59 

0.33 

0.59 

0.44 

0.31 

0.52 

0.29 

PL-ECT3 

0.58 

0.54 

0.61 

0.54 

0.27 

0.17 

0.30 

0.18 

TG-ECT 

0.27 

0.31 

0.23 

0.31 

0.67 

0.67 

0.67 

0.67 

Table  6 
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